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' Abstract The Backlund transformation for pseudospherical surfaces, which is equiv- 

^\ , alent to that of the sine-Gordon equation, can be restricted to give a transformation 

on space curves that preserves constant torsion. We study its effects on closed curves 
(in particular, elastic rods) that generate multiphase solutions for the vortex fila- 
^ ' ment flow (also known as the Localized Induction Equation). In doing so, we obtain 

, analytic constant-torsion representatives for a large number of knot types. 
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Introduction 

Soliton equations have become a familiar presence in the differential geometry 
00 ! of curves and surfaces. The description of pseudo-spherical surfaces and their 

^ I asymptotic lines in terms of the sine-Gordon equation dates back nearly a century. 

Q\ I Of roughly the same date is the derivation (see [Ri]) of the Localized Induction 

' Equation (LIE) , which provides one of the richest examples of connection between 

curve geometry and integrability. 
bp! The understanding of this connection has progressed in recent years along dif- 

ferent directions. On the one hand, several fundamental properties of soliton 
equations have been given a geometrical realization; in the case of the LIE, its bi- 
^ ■ hamiltonian structure and recursion operator, its hierarchy of constants of motion, 

and its relation to the nonlinear Schrodinger equation possess a natural geometric 
interpretation [L-Pl]. 

On the other hand, some well-known classes of curves in differential geometry 
have been identified with solutions of integrable equations: for example, elastic 
curves and center-lines of elastic rods are among the solitons for the LIE, curves of 
constant torsion correspond to characteristics for the sine-Gordon equation, and 
planar and spherical curves are associated with solutions of the mKdV hierarchy 
([G-P],[L-P2],[D-S]). 

A third direction of research ([C],[M-R]) concerns the topological properties of 
closed curves that arise as solutions to soliton equations. A major question is 
whether the presence of infinitely many symmetries and the associated sequence 
of integral invariants may be related to knot invariants; in fact, one hopes that 
the knot types of the periodic analogues of soliton solutions (the so-called multi- 
phase solutions) can be described using methods of integrable systems, such as 
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the Floquet spectrum and Backlund transformations. We remark that the knot 
types of closed elastic curves, which turn out to generate two-phase solutions of 
the LIE, were classified by Langer and Singer [L-Sl] with the aid of infinitesimal 
symmetries related to flows in the LIE hierarchy. 

In this work we illustrate these different approaches through the study of a uni- 
fying example. First, we propose a geometric realization of the Backlund transfor- 
mation for the sine-Gordon equation in the context of curves of constant torsion, 
by restricting the classical Backlund transformation for pseudo-spherical surfaces 
[C-T] to a constant-torsion-preserving transformation between asymptotic lines. 
Because it is not necessary to embed the curve in a surface to carry out the 
transformation, this is also a tool for obtaining new interesting curves of constant 
torsion. 

Secondly, we are interested in curves that generate special solutions (in partic- 
ular, solitons and their periodic counterparts) of integrable equations and whose 
analytic expressions and general properties are well-known. This choice is par- 
ticularly significant for closed curves, since multi-phase solutions are dense in the 
space of all periodic solutions: it is therefore important to understand their ge- 
ometrical and topological properties. The simplest non-planar constant torsion 
curves in this class are the centerlines of elastic rods. (Elastic rods, which also 
generate two-phase solutions of the LIE, are critical curves for the elastic energy 
(total squared curvature), subject to fixed total torsion and fixed length.) 

Finally, we study the topological properties of curves obtained by means of 
(single and iterated) Backlund transformations of constant torsion elastic rods. 
The behaviour of these transformations is diverse: we show that a single Backlund 
transformation preserving closure leaves the class of constant torsion elastic rods 
invariant, and we compute an interesting formula relating the linking number of 
a curve with its transform and the self-linking of the original curve. The iter- 
ated Backlund transformation instead allows one to leave the class of elastic rods 
and to produce a variety of phenomena such as knotting, self-intersections and 
unknotting. Effects of Backlund transformations on the multi-phase solutions of 
the sine-Gordon equation have been discussed in the work by Ercolani, Forest and 
McLaughlin [E-F-M] , which uses a combination of Floquet theory for the associ- 
ated spectral problem and methods of algebraic geometry. We present here the 
first concrete realization of their conclusions in a geometric setting, and give simple 
direct proofs for the case of elastic rods of constant torsion. 

This article is structured as follows. In order to implement the Backlund trans- 
formation for pseudospherical surfaces restricted to constant torsion curves, we 
need to solve an associated Riccati equation. To overcome this difficulty, we prove 
in Part I that, if we know the expression of the Prenet frame of a given curve that 
depends analytically on the constant torsion r, then we can obtain the solution to 
the Riccati equation by analytic continuation to an imaginary value of the torsion. 
In other words, for any curve whose Frenet frame can be continued analytically 
as a function of r, with fixed curvature function, the Backlund transformation 
is directly computable. This is illustrated for closed elastic rods of constant tor- 
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sion. Using Langer and Singer's solution of the Euler-Lagrange equations [L-S4], 
we derive explicit formulas for these curves and their Frenet frames, and we also 
describe their knot types. Wc then show that the condition that the Backlund 
transformation give another closed curve forces the new curve to be an elastic rod 
that is congruent to the original. However, such a closed Backlund transformation 
carries topological information by providing a measure of the knottedness of the 
initial elastic rod. In §1.6 we show this by relating the linking number of the curve 
with its transform to the self-linking of the original curve. Part I ends with a 
discussion of, and a precise conjecture about, the effect of Backlund transforma- 
tions on constant torsion n-phase (and n-soliton) solutions of the LIE hierarchy. 
The conjecture, which is easily verified for the first few values of n, concerns the 
form of certain Killing fields associated with n-soliton curves, and it implies that 
n-solitons are in general taken to (n + l)-solitons. 

Part II concerns the iterated Backlund transformation, which consists of two 
successive Backlund transformations using related solutions of the Riccati equa- 
tions. For the purpose of constructing iterations, we conveniently rederive the 
original Backlund formula in terms of a gauge transformation of the sine-Gordon 
linear system. (This is simply the spatial part, in characteristic coordinates, of 
the Lax pair, and is equivalent to our Riccati equation.) Using this formulation, a 
simple algebraic procedure produces the iterated formula which we apply to con- 
stant torsion elastic rods to produce a variety of interesting curves. The closure 
condition, developed using standard arguments of Floquet theory, is in this case 
less restrictive. As a consequence, we prove that, beyond producing congruent 
elastic rods, an iterated Backlund transformation can take constant torsion two- 
phase solutions to multi-solitons whose knot type differs dramatically from that 
of the original curve. All three phenomena occur: knotting, self-intersections and 
unknotting, all related to the Floquet spectrum of the linear system associated to 
the initial curve, as illustrated and discussed in the concluding section. 

The investigations that we report here were initiated by, and advanced with the 
steady encouragement of, the geometry group at Case Western Reserve University 
(Joel Langer, David Singer, and the authors), with important contributions from 
Ron Perline of Drexel University. In particular, Perline first asked if there exist 
constant torsion solitons for the LIE; Singer classified the elastic rods of constant 
torsion; and, Langer and Perline developed the theory of planar-like solitons and 
their associated Killing fields which allows us to connect our Backlund transfor- 
mations with the LIE hierarchy. We gratefully acknowledge these contributions, 
and the lively discussions that helped to shape this work. 

I. The SINGLE Backlund transformation 

1.1 Backlund transformations and the Frenet frame 

In classical differential geometry, a Backlund transformation takes a given pseu- 
dospherical (i.e. constant negative Gauss curvature) surface to a new pseudospher- 
ical surface. As explained by Chern and Terng [C-T] , the new surface is connected 
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to the old surface by line segments that are tangent to both surfaces, of a fixed 
length, and such that the angle between the surface normals at corresponding 
points is also constant. Moreover, the Backlund transformation takes asymptotic 
lines to asymptotic lines. Since the asymptotic lines on a pseudospherical surface 
have constant torsion, it is not surprising that we can restrict the Backlund trans- 
formation to get a transformation that carries constant torsion curves to constant 
torsion curves. 

Theorem 1.1. Let j{s) be a smooth curve of constant torsion r in ^ , parametrized 
by arclength s. Let T,N,B be a Frenet frame, and k{s) the curvature ofj. For 
any constant C , let j3 = (3{s; k{s), C) be a solution of the differential equation 

^ = CsmP-K. (1) 
ds 

Then the curve 

2C 

lis) = lis) + -^^-^ioosPT + sin (3N) (2) 
is a curve of constant torsion r, also parametrized by arclength. 



Proof. Since the angle between the normals of the two pseudospherical surfaces is 
constant at corresponding points, it follows that the angle between the binormals 
of 7 and 7 is constant. In fact, one computes the Frenet frame along 7: 

f = T + (1 - cos ^) sin /3(cos (3N - sin f3T) + sin 9 sin (3B 

iV = - (1 - cos ^) cos /3(cos I3N - sin /3T) - sin 9 cos I3B 
B = cos dB + sin 6'(cos /3iV - sin /JT), 

where tan(^/2) = C jr. Then the torsion is r and the curvature is 

/5 = Av-2Csin/?. (3) 

□ 

We should point out that we will assume our curves have a smooth generalized 
Frenet frame, that is, smooth orthonormal vector fields T, N, B that satisfy the 
Frenet equations with respect to smooth functions k, t. In particular, k is allowed 
to change sign. Given an orientation of 7, T is well defined, and A^, B, k are well- 
defined along such a curve up to multiplying all three by —1. If the curve is closed 
of length L, then our assumptions imply that there must either be a smooth Frenet 
framing by vector fields of period L, or a framing by vector fields of period 2L 
such that N,B, K change by a minus sign after one circuit of the curve. (We will 
refer to these closed curves as 'even' and 'odd', respectively.) 



Example. To compute the Backlund transformation of a constant torsion curve, 
one needs to solve the differential equation (1), which is equivalent to a Riccati 
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equation when one changes variables to y = tan(/3/2). When k is equal to a 
constant kq, the curve 7 is a helix, and the ODE is easy to solve. For C < kq, 
the solutions are periodic, and one can arrange that the period is commensurable 
with the translational period of the helix, yielding a curve 7 that is periodic up 
to translation (see Figure 1.1(b)). For C > kq, the solutions are asymptotically 
constant, and yield curves 7 that are asymptotic to a helix (see Figure 1.1(c)). For 
C = kq, the curvature of 7 is a rational function of arclength: 




Figure 1.1. A helix and its periodic and asymptotically helical Backlund 
transformations 

The Backlund transformation for pseudospherical surfaces is, of course, equiv- 
alent to the Backlund transformation for the sine-Gordon equation. If s and t 
are arclength coordinates along the asymptotic lines, and 9 is the angle between 
the asymptotic lines, then the Codazzi equations imply 9xt = sin^. In terms of 
(p = 9/2, the Backlund transformation for sine- Gordon is 



{(fi-(p)8 = Csin{(p + (p), {(fi + (p)t = C ^sm{(p-(p). 
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If T, N are chosen along an s-curve so as to agree with the orientation d/dt Ad/ds 
on the surface, then ds = k(s). Once we make this identification, then (1) follows 

Along with applying the Backlund transformation to constant torsion curves, 
we can also apply the nonlinear superposition principle for solutions of the sine- 
Gordon equation [Ro]. (For surfaces, this is Bianchi's theorem of permutability 
[E].) Suppose (3i = /3(s; k(s), Ci) and P2 = «;(s), C2) generate two Backlund 
transformations, with corresponding curvature functions Ki{s), ^2(5) given by (3). 
Then particular solutions (3i2 = C2) and P21 = P{s', K2{s),Ci) to suc- 

cessive Backlund transformations are given by 



/3i2 - /3i = /321 - /32 = 2 arctan 



C1 + C2 
Ci — C2 



tan 



/3i-/32 



(4) 



Now suppose that we can solve (1) for some initial curvature function and all 
values of C, as we can for the helix. Then in principle we can use (4) to calculate 
arbitrarily many Backlund transformations of the initial curve, choosing constants 
Ci, C2, C3, . . . . (Recall also that once a single solution of a Riccati equation is 
known, all others can be obtained by quadrature.) The superposition principle 
(4) also makes sense in the limit as C2 ^ Ci, if we fix initial values for /?i and P2 
independent of Ci, C2. 



1.2 The Frenet frame and the linear system 

In this section we will see that, if we know the solution of the Frenet equations 
for a given analytic curvature function k{s), as analytic functions of the constant 
torsion r, then we already know the solution to the Riccati equation (1) that 
we have to solve to compute the Backlund transformation for the corresponding 
constant torsion curves. 

For the given curvature function k(s), suppose '^{s; A) is the fundamental solu- 
tion matrix for the linear system 



(^_lf A k{s) 
ds ~ 2 V -At(s) -A 



V'. (5) 



(It follows from standard ODE theory that \E' is an analytic function of A.) When 
A = —IT for T e M, this system encodes the Frenet equations, in the following way. 
If 

/ N -T-l 



then 7 is a curve in stt(2), which we will identify with (see below). One can 
verify that 7 satisfies the Frenet equations for curvature k(s) and constant torsion 
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T, when the Prenet frame is given by the su{2) matrices 



2 \0 -i 
^=i*-^f- (6) 



2 \i 

Now suppose that ^'(s; A) were extended complex analytically from A pure imag- 
inary to all A G C. We observe then that for A = C G M, the coefficient matrix 
in (5) is real, and for any real vector solution i/j, the ratio of entries y = — 'i/'i/V'2 
satisfies the Riccati equation 

^ = C^2/- 2(1 + 2/). 

Since this is the same as the equation satisfied by tan(/9/2) under the Backlund 
transformation, we obtain a Backlund transformation for the curve 7 using 

In practice, we may not know the fundamental matrix \1/ explicitly, but if we 
know the Frenet frame for the curve, normalized so that (T, N, B) is the identity 
matrix at s = 0, and use the identification 

(i,o.o)^l(; _°^). (o.i.o)^l(° (o.o.i)^l(j 

of with sw(2), then (6) gives us values for all possible quadratics in the entries 
of \E'. Thus, for any linear combination ij) of the columns of \E', we can express the 
quadratics in (7) in terms of analytic extensions of those obtained using (6). 

1.3 Torus knots and closed elastic rods of constant torsion 

In a recent paper [L-S4], Langer and Singer formulate and solve the Euler- 
Lagrange equations associated to the Kirchojf elastic rod, by which is meant a 
curve that is critical for some linear combination of J n'^ds, J rds and length. 
Among the solutions they obtain are two families of constant torsion curves. For 
the family we will use, the curvature is given by an elliptic cosine, and the shape 
of the rod is governed by the torsion, the maximum curvature kq, and the elliptic 
modulus p: 

k{s) = Kocn {x,p), where x = 
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(In what follows, we will take kq = 1.) Up to scale, then, the shape of the rod 
depends only on two parameters; we will use the parameters p and a = 2t, and 
the coordinate x along the rod. We will show how to choose the parameters to 
obtain closed elastic rods of constant torsion. 

Langer and Singer compute the position of the rod in a system of cylindrical 
coordinates (r, 0, z) that are generated by the Killing fields associated to the rod 
(see [L-S4], §5). In these coordinates, r is already a 2iC-periodic function of x, 
where K = K{p) is the complete elliptic integral of the first kind: 



r = —\ sn^a;, 

/X V TO 

where 

m = IGfi^/ip-^ + cr'^f. 



Closure in z imposes (cf. equation (26) in [L-S4]) 



.2 1 ,„E{v) 



where E{p) is the complete elliptic integral of the second kind. (David Singer has 
shown that constant torsion rods in the other family, whose curvature is given 
elliptic dn, are never closed in the z coordinate.) This implies |cr| < p~^ and that 
p has a maximum possible value Pmax ~ .9089085 at which approaches zero. 
(We will take a > 0, giving positive torsion curves; curves obtained using a < 
will differ by a refiection.) Using the remaining parameter p, we can attempt to 
make the change in 9 over a 2K period equal to a rational multiple of 27r. The 
Euler-Lagrange equations imply that 

dO per (p-2-l)(p-2-(T2) \ , . .2 2 .N/. 

— = -— Ai + f- , where Ai = (a^ - p'^ + 2 /4. 

ax \ Z{p ^ + a"^j(l — msn ^x) / 

(This formula is due to Singer.) So, the change in 6 can be expressed in terms of 
complete elliptic integrals: 

A6' (E \ Xipa 



2K 



where 
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(cf. formula 413.01 in [B-F]). 

Remark. It will be desirable to have a formula for A(p, a) that does not involve ^ as 
denominator. Using Legendre's relation, and the addition formula and imaginary 
transformation for Jacobi's zeta function Z{x) (see [W-W], §22.73), we obtain 

A = —iZ (sn~"'^Q;) +p^crQ!, where a = — — tJ^ — ^. (10) 



The following result now follows directly from the above formulas and those in 
[L-S4]: 

Theorem 1.2. Given any relatively prime integers m,n such that \m/n\ < 1/2, 
there exists a smooth closed elastic rod of constant torsion with the knot type of an 
(m, n) torus knot. 

Proof. We use continuity, noting that the limiting values of A^/27r are zero as 
p\0, and one half as p — > Pmax, to assert that A^/27r = m/n for some modulus 
p. Then one checks that, over the course of one 2K period, r{x) and z{x) describe 
a simple closed curve in the rz plane. □ 

Remarks. It is interesting that this is exactly the same set of knot types as 
are available among elastic curves (see [L-Sl]), which of course are critical for a 
combination of J K^ds and length, and do not have constant torsion (unless they 
are planar). Since a computer-generated plot of shows it to be a monotone 
increasing function of p, we expect that, similar to the situation for elastic curves, 
there is a unique constant torsion elastic rod for each possible torus knot type. 
(However, by arranging that A0/27r = 1/n for any n > 3, we can get many 
unknots.) Some of these torus knots, as well as some unknots, are shown in Figure 
1.2. It is also apparent that, as p — > Pmax and the torsion approaches zero, the 
curve approaches a figure-eight elastic curve contained in a plane through the 
2;- axis. 



For the purpose of computing their Backlund transformations, we need to com- 
pute Frenet frames for these elastic rods. Using the formulas in [L-S4], one can 
express the Frenet frame in a cylindrical coordinate basis: 



T 

N 



( cn,Tsn,Tdna; 

cna;(p~^ -\- — 2sn^a;) 
(jsn xdn x 



g^3^2 

Aicjsn xdnx 



2pii 



a{p ^ -I- - 4Aisn ^x) cn^a;-2Ai\ 

2/1 
sn a;dn x 

2pii 
acnx 

2^ J 



2pn^r'^ 



16fi^r 



a )cnx 



To convert to Cartesian coordinates, we need the cylindrical coordinates r, 6 as 
functions along the elastic rod. Using formulas that express incomplete elliptic 
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Figure 1.2. Some closed elastic rods of constant torsion. (The curves 
here have been thickened by a vertical ribbon.) 

integrals of the third kind in terms of theta functions (cf. formula 434.01 in [B- 

F]), we get 



r cos 6 



rsYn.9 



2Kp' e-'^"" Qiix - iF) + e'^^'Qiix + iF) 
Trp 2^Q{x)Hi{iF) 



2Kp' e-*^^ei(x-iF) -e*^^ei(a; + zF) 
T^P 2fiie{x)Hi{iF) 



where F = F{^,p') and A = A{a,p) is defined in (9). 

Now let H be the matrix whose rows are T,N,B. Let Hq = H\j:^o, and let 
G{x,a) = Hq^H; G is the normalized Frenet frame. The above formulas can be 
used to express the entries of G as complex analytic functions of a in a disk about 
the origin in the cr-plane. (One can also desingularize the formulas near the points 
(T = ±i ± p' /p, where n — 0.) In particular, from G{x, 2iC) we may, using the 
formulas (7) in the previous section, construct the Backlund transformation of an 
elastic rod of constant torsion. However, the new curve is not necessarily closed; 
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in subsequent sections we will see how to obtain closure and how severely closure 
restricts the shape of the new curve. 

Remarks. 

1. From the formulas in [L-S4] one can also obtain 

= ^§4. (12) 

This b striking resemblance to the formula obtained by Mumford [M] for 

planar elastic curves in terms of theta functions. In fact, one can check that as 
P Pmax and our elastic rod becomes planar, (11) and (12) agree with Mumford's 
formula in the limit. 

2. One can derive formulas quite similar to (11) and (12) for elastic rods in 
general. 

1.4 Closure of the Backlund transformation 

Suppose we have generated a Backlund transformation of a constant torsion 
curve 7, by obtaining the fundamental matrix ^ for the linear system (5), selecting 
a non-zero vector solution 

V' = *(^j) ci,c2eM, 

and performing the transformation using (7) and (2). Suppose that 7 is closed 
and has length L. Then the Backlund transform 7 will be closed if the vector 
V = cos (3T + sin (3N has period kL for k a positive integer. 

Proposition 1.3. Suppose V has period kL. If k is even or 7 is even, then 
(ci,C2) is an eigenvector o/^'Is^fcL- If k is odd and 7 is odd, then (ci,C2) is an 



eigenvector of (^^ *|s=fcL- 



When 7 is an elastic rod of constant torsion, we may use the entries of the 
normalized Frenet frame G{x,a), together with the formulae (6), to obtain all 
quadratics in the entries of ^. In particular, to obtain the eigenvectors required 
in the above proposition, we may use the squares of the relevant matrices. 

Suppose that the elastic rod closes up after n periods of length 2K in the x 
parameter. (The rod is even or odd exactly as n is even or odd.) If kn is even, 
then we calculate that 



^'^\x=2knK= cos{xA) - sin(a;A) 



2// 2/x 
_a . 1 - 2Ai 

V ~2/I ~* 2/1 ) 



where A is given by (9). Note that, when we want to calculate Backlund transfor- 
mations, we set (7 = 2t = 2iC for C e R, C 7^ 0, and then A is purely imaginary 
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and is real. Moreover, the eigenvalues cos(a;A) ± i sin(a;A) of this matrix are 
always distinct. If kn is odd, then 



( ( -1 ) *) \x=2knK= - cos(a;A) ^ J J ^ + sin(a;A) 



/.1-2A1 o_ X 



2/x 
cr 

V 



2/1 

1 - 2Ai 

2iJ, ) 



— i.e., the same matrix as above, up to a minus sign. Thus, in either case the 
eigenvalues are distinct. The corresponding eigenvectors are 



v+ = 



V 



2/i 

1 -2Ai 



2/1 ) ) 



V- = 



V 



a 
2/1 



(13) 



Since these eigenvectors are independent of k, we conclude that when is a dosed 
elastic rod of constant torsion, 7 closes up after k circuits around 7 if and only if 
it closes up after one circuit. 

1.5 Elastic rods again 

In this section we will see that, for any value of the Backlund parameter C, a 
closed Backlund transform of a closed elastic rod of constant torsion is congruent 
to that rod. Since the two curves have the same torsion, it will suffice to show 
that they have the same curvature function, up to a phase shift. 

Let q = 2C and a = iq ior q e IR. Then by (7) the initial values for P 
corresponding to eigenvectors v+ and V- given above are 



sin/9 



2q 



COS/? = ±^^^ . 



p-2 _|_ q2 ' 

Since K = cnx, the differential equation for /? is 

d(3/dx — p{q sin (3 — 2cnx) 
If K were cn (x — a) then (3) would imply that 

q sin P = cnx — cn{x — a) 

and 



(14) 



(15) 



(16) 



dj3/dx = —p{cnx + cn (a; — a)). 



Then by integration. 



/? = — [sin ""^ (p sn x) + sin ^(psn(x — a))] + Ci 



(17) 
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for some constant Ci , giving 



sin(/? — Ci) = —p[snxdn (x — a) + dnxsn (x — a)] 
cos(/3 — Ci) = dna;dn (x — a) — p^sna;sn (x — a). 

We can find the constants a and Ci by setting a; = 0: (16) gives 

cna = {p~^ - q^)/{p~'^ + q^) 

and so 



2p 
p-2 _|_ q2 

while (17) gives sin(/3 — Ci) = psna when x = 0. Comparing with (14) shows that 
we should use Ci = or Ci = tt; since cos(/3 — Ci) is always positive in (17), we 
use Ci = and the plus sign in (18) for v+, Ci — n and the minus sign in (18) for 

V-. 

It remains only to be seen that (17) gives a solution to (15). (We will show this 
for the initial values associated to using Ci ~ and the plus sign in (18); the 
V- case is similar.) Let b = a/2. Then from (17), 



dp/dx = —2pcnx — 2p 
(cf. 123.03 in [B-F]). Since 

pq s'mP = —2p"^q 



sn 6 dn 6 sn (x — b)dn. {x — b) 
1 — p'^sn ^6 sn ^ (a; — b) 



2 sn (a; — 6)dn (x — b)cn b 
1 — p'^sn^bsn^{x — b) ' 



we will be done if we can show that pqcn b = sn bdn b. This now follows from the 
double angle formula 

sn a sn 6 dn 6 



1 + cn a cn 6 

Now, by construction, k = cn (a; — a) and hence 7 is congruent to the elastic rod 
7- 



1.6 Linking numbers 

In this section we use the explicit solution of (1) obtained for an elastic rod 7 in 
the previous section, to compute the linking number of 7 and its twin 7 in terms 
of the self-linking of 7 (see [Pol] for discussion of self- linking) . 

Theorem 1.4. If \C\ is sufficiently small, 



Lk{^,^) = SL{^)-n/2 



(19) 
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where n is the number of 2K periods of ^. 

The proof will be an application of White's formula (see [W], [Po2]). Before 
proving this result, we should remind the reader that the self-linking number 
SL{'y) can be thought of as the linking number of 7 with the curve traced out 
by the endpoint of the Frenet normal (assuming 7 is even), suitably scaled so 
that the vector never intersects other points of 7. As such, SL{'y) is the sum 
of a contribution from the twisting of N about 7 and a contribution from the 
twisting of 7 about itself — a measure of the knottedness of 7. In (19), the former 
contribution is cancelled out by the —n/2. In particular, when 7 is unknotted (as 
it is when, for example, we arrange that = 7r/2), then LA;(7,7) = 0. 

Proof. Let us recall White's Formula. Let 7 be a smooth closed space curve, 
oriented by unit tangent vector T, and V a unit normal vector along 7. Let 7^ be 

the curve traced out by the endpoints of SV, with the obvious orientation, where 
5 > is chosen small enough that the ribbon spanned by SV, with boundary 
7 [J 7^, is embedded. Then 

Lfc(7, 7^) = I^r(7) + ^J{TxV)-dV, (20) 

where the writhe is given by the integral 

Wr{^) = j j e*dS, 
7x7 

wherein for (x, |/) G 7 x 7, a; 7^ y, e is the unit vector from x to y, and dS is the 
element of area on the unit sphere. Before White, Pohl [Pol] obtained a special 
case of this formula, 

SL{^) = Wr{^) + Tds, (21) 

which allows us to express the writhe in terms of the self- linking of 7. (While Pohl 
did not define the self- linking number for curves with infiection points, we will 
take (21) as defining SL{'y) for such curves. If 7 is an even curve, this SL{'y) is 
indeed the linking number of 7 and 7jv — for either choice of N — while if 7 is odd, 
SL{'y) is a half-integer.) 

We next remark that it is easy to extend (20) to the case where F is a transverse 
vector field along 7: 

LA;(7, 7^) = Wr{^) + 1- J (T x U) ■ dU, (22) 

where U is the unit vector in the direction of the orthogonal projection of V into 
the normal plane along 7. However, for the Backlund transformation we have 



F = cos /? T + sin /? A^, 



BACKLUND TRANSFORMATIONS AND KNOTS OF CONSTANT TORSION 15 



which fails to be transverse to 7 whenever sin /5 = 0. In order to use (22) to calcu- 
late the linking number, we will take a small perturbation VofV that is transverse 
to 7, also ensuring that the ribbon spanned by V is embedded. Then, since, as 
shown by the Gauss integral formula, the linking number depends continously on 
the two curves, our calculation of L/c(7,7^) will give L/c(7,7^) also. 

We may assume that k and r 7^ in the vicinity of any point where V is 
tangent to 7. A local calculation then shows that, when F is a positive multiple 
of T at the point of tangency, the perturbation 

V = V - €ktB, € > (23) 

makes the ribbon spanned by dV embedded, for d and e sufficiently small. (It 
is interesting to note that if e is negative in (23), then no matter how small 5 is 
chosen, the ribbon fails to be embedded near the point of tangency. The reader 
may visualize this by imagining a momentarily tangent vector field V along a 
right-handed helix, with k > and r > 0: the vector field must be perturbed 
down, not up, in order to prevent the perturbed vector field from intersecting the 
helix.) Now (22) gives 



LA;(7,7 ) = SL{'y) + 77- / rfarctan 
V 2n J 



sin/3 



eKT 



When 7 is a constant torsion elastic rod, and the Backlund transform is con- 
structed using eigenvector v+, 

. ^ cnx — cn(x — a) 2pq 

sm/J= ^ -, sna=- ^ 

q 1 -l-p^g^ 

and cos/3 > always, showing that our assumptions hold. Furthermore, sin/3/K 
goes from -|-oo to —00 between infiection points, and hence 

LA;(7,7y) = '5L(7)-n/2. 

The same result holds when we use eigenvector v_ , but since cos /? < we have to 

use e < in (23); however, sna = ^-^ makes up for the change in sign. □ 

1 + p^q^ 

1.7 Backlund transformations and the LIE hierarchy 

The constant torsion elastic rods are among the "soliton" curves for the Local- 
ized Induction Equation (LIE), and as such arc critical for some linear combination 
of the integral invariants associated with the LIE. In this section, we will formulate 
a precise conjecture as to how the Backlund tranformation takes solitons to soli- 
tons, and in particular how the coefficients in the linear combination of invariants 
change. Ercolani, Forest and Mclaughlin [E-F-M] use arguments from algebraic 
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geometry to address the analogous question for n-soliton solutions of the sine- 
Gordon equation. They find that in general a Backlund transformation (3) takes 
n-solitons to (n-l-l)-solitons; what we describe below is a more geometric argument 
supporting the same conclusion. In the particular case when the initial condition 
of the Riccati equation selects one eigenvector of the transfer matrix of the linear 
system, producing closed curves as we discussed above, then [E-F-M] find that 
the Backlund formula takes an n-soliton to another n-soliton; this of course agrees 
with our results on closed Backlund transformations of elastic rods. 

We begin with the observation, due to Langer and Perline [Pe], that along a 
curve of constant torsion r the linear combinations of the LIE vector fields {X^j'^Q 
defined by 




n > 



are purely binormal and do not involve r. For example, 

Bo = Xi = kB, Bi = X3- 2tX2 + T^Xi = -{k" + k^/2)B. 

A planar-like n-soliton [Pe] is a curve of constant torsion along which some constant- 
coefficient linear combination ^^^q '^kBk vanishes. For example, a helix is a 
planar-like 1-soliton. Because the LIE vector fields X^ are Hamiltonian vector 
fields associated with certain densities involving k, t and their derivatives (see 
[L-Pl] for details), these solitons are critical curves for a natural variational prob- 
lem. For example, since Xi, X2, X^ are Hamiltonian vector fields for J ds, J rds, 
J ^k^ds, an elastic rod of constant torsion is also a planar-like 1-soliton. 

It is easy to show, using the LIE recursion operator [L-P2], that if along a 
constant torsion curve 7 we define 

2n 
k=0 



2n 
k 



i-r) 



2n — k V 



(24) 



then X^^^o (^kBk = if and only if 



n 

J = J2 ""kUk (25) 

k=0 

is a constant vector. We call J a Killing field along 7 because it is the restriction of 
a vector field generating a symmetry of Euclidean space, in this case translation. 
Combinations of LIE vector fields that are Killing fields play an important role in 
showing that variational problems associated to combinations of the LIE function- 
als yield finite-dimensional completely integrable Hamiltonian systems (see [L-S2] , 
[L-S3], [L-S4]). So, it is important to see how Killing fields along planar-like n- 
solitons interact with the Backlund transformation. 
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Conjecture 1.5. Let ^ be a constant torsion n-soliton, and let ^ be a Bdcklund 
transformation of^. Let Uk denote the combinations of LIE vectors defined by 
(24), computed along 7. Then if J is the constant vector defined by (25) along 7, 

n n 
.k=0 k=0 

where C is the parameter in (1), (2) and h is a constant depending only on 7. 

For any given n, one can mechanically verify this conjecture using the formulae 
for the vector fields Uk-, Uk-, the equation (3) for the curvature of 7, the differential 
equation (1), and the ODE for the curvature of 7 implied by Y^=Q^kBk = 0. 
(The constant h arises as a first integral of this equation.) 

We have verified the conjecture forn = 0, 1, 2; hence we know that the Backlund 
transform takes a planar-like 1-soliton to a planar-like 2-soliton, and takes a planar- 
like 2-soliton to a planar-like 3-soliton. However, the results of previous sections 
show that a closed transform of a closed 1-soliton is still a 1-soliton; an intuitive 
explanation for this can be given for small values of the arbitrary constant C in 
the Backlund formula. When C — 0, since sin (3 is bounded, the Riccati equation 

dl3 

becomes approximately — ~ —n. Therefore (3{s) ~ — ^'(-s) = — K{s')ds' and, 
to first order, 

2(7 

7 = 7 + -^ [cos er - sin ON] . 

In other words, the Backlund transformation flows along the vector field 

7^ ~ = cos 6^ - sin 6'iV. 

This vector field W, discussed in [L-P2], is known as the trigonometric vector field, 
and the associated evolution for = K{s')ds' is described by the sine-Gordon 
equation, which preserves soliton type. Moreover, for periodic solutions of the 
sine-Gordon equation, the limiting Backlund formula described above transforms 
closed n-soliton curves into closed n-soliton curves. 

In the second part of this paper, we will perform iterated Backlund transfor- 
mations on the 1-solitons to obtain new closed curves that are multi-solitons for 
the LIE hierarchy. 

II. The double Backlund transformation 

In Part I we used Backlund transformations to produce new constant torsion 
elastic rods from given ones. Moreover, we verified that imposing the closure 
condition on the transformed curves restricts them to the same congruence class 
as the original elastic rod. 

More exotic and interesting constant torsion curves are produced by succes- 
sive iterations of the Backlund transformation starting with a given elastic rod of 
constant torsion. 



t2 + C2 
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2.1 A GAUGE FORMULA FOR THE BACKLUND TRANSFORMATION 

We present an alternative derivation of the Backlund formula based on a gauge 
transformation for the associated linear system, as described in [E-F-M]. This 
approach is most convenient for constructing iterated Backlund transformations, 
since it reduces each step of the iteration to a purely algebraic computation. We 
will show that the Riccati equation ( equivalent ly, the linear system) needs to be 
solved only once at the initial step of the iteration. 

We begin with the following form of the associated linear system 



ds 



1 / iK{s) 

2 V A 



A 

—iK{s) 



(26) 



which simplifies the formulation of the Backlund transformation. We observe that 
equation (26) is transformed into the linear system (5) by setting = Aip, where 
1 

Using the fundamental solution matrix $ of (26), we identify the curve of con- 
stant torsion r with the su{2) matrix 



7(s) = 

and compute its Frenet frame 



.1 

dX 



X=—iT 



(27) 



T 
N 
B 



1 
2 
1 



2 



i 

1 





1 




i 



(28) 



The Backlund formula is a consequence of the following 

Proposition 2.1. Let x(s; i/) be a solution of the linear system (26) for X = v, 
and let </>(«; A) be a solution of (26) for an arbitrary A e C. Then, the gauge 
transformation 

A -V- 



-V- 



,X2 

Xi 



produces a solution (f)^^^ of 



X2 



X 



</>(«; A), 



(29) 



dcf) 
ds 



1 
2 



■mW(s) 



A 

-mW(s) 



0, 



(30) 
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where 

(31) 



Vx2 xij 



The proof is a direct verification of formulas (29) and (31). 

We will note for future reference that z = Xi/X2 satisfies the Riccati equation 
dz/ds = iKZ + v{l - z'^)/2. 

The gauge transformation (29) can be normalized to produce the fundamental 
matrix for (30): 

Assume that k{s) is real, and suppose for the moment that A = —ir for r e M. 
The transformation (31) does not, in general, produce a real-valued function k*^^-* (s) 
unless some condition is imposed on v. Note that the matrix $ now takes values 
in SU{2), and (32) becomes 

This gauge transformation is in SU{2) only if uz = v/z, which implies that v is 
real and z has unit modulus. In this case, k*^^^ is real. Using the reconstruction 
formula 

7(s) = ' 



dX 

we obtain the following expression for the transformed curve: 

If we make the substitution z = —e~^^, then (3{s) satisfies the equation (1) for 
the Backlund transformation, with v = C, and (33) becomes (2), in terms of the 
Frenet frame for 7 given by (28). 

2.2 The iterated formula 

The gauge transformation (29) can be iterated algebraically once the solution 
of the linear system (26) at an initial curvature function «:(s) is known. In fact, 
formula (29) produces the eigenfunction of (30) out of which the new gauge matrix 
is constructed. 

We begin with some symmetry considerations. Given a solution x(s, v) of the 
linear system (26) at a pair (k(s), v)-, then 

X(^,^)=(J j)xM. (34) 
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is a solution of (26) at v). We use this property to describe the following 

iterated Backlund transformation 



7 



7 



(1) 



7 



(2) 



(The use of the symmetry (34) is crucial in constructing real-valued curvature 
functions and curves in su{2).) The following lemma is obtained by direct compu- 
tation: 

Lemma 2.2. Let f{s; i>) be a solution of (26) at {K{s),i>), and let 

Xi 



(is; 



( 



V 



-V- 



X2 

t 

XI 



X2 



he the new solution at {k^{s), u) constructed by means of the gauge transformation 
(29) . Then, two applications of the gauge formula (29), using solutions x(s; zv) 
and C,{s;V), gives 



1 



X1C2 X2 C2 

xi Ci X2C1 / 



>;A). 
(35) 



Moreover, the matrix T{—iT, v), t eM., is unitary if and only if we choose f{s,v) = 

(J j)xM- 

We can now derive the explicit formula for the transformed curve and for its 
corresponding curvature function by using the reconstruction formula (27) and 
the form (35) of the new eigenfunction. We set Xi/X2 = pe~'^ (no longer unit 

modulus), A = —IT and we introduce the constant a = i—^ ^, , „ 

(r^ _l- zy^j(^r^ -I- z/^J 

Proposition 2.3. The iterated Backlund transformation of the curve 7(5) at the 
eigenvalues v and v is the constant torsion curve 

7(2) (s) = ^(s) + « (pT + ^AT + rB) , 



(36) 



where 
V = 
9 = 



-1 



{y — vp^){y — vp^) 
1 

{v — up'^){u — z/p2) 



V — V 

2i 



(r^ + + p2)pcos/5+ ^(r^ - - p^)psmp 



^(r2-M^)(l-p>cos/3 



2 

f — u 



2i 



— {t'+\u\^){1 + p^)p sin P 



r = Tiui 



— vp'^){v — vp"^) 
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The associated curvature function is 

K^'Hs) = k{s)+ . ^^''i';' [iu + P)(p2 - 1) COS/? + z(P - u)ip^ + 1) Sin/3] . 

(37) 

Remarks. 

1. We observe that the expression for 7*^^) — 7 now contains a binormal compo- 
nent proportional to the torsion of the original curve. 

2. The double Backlund transformation makes sense for all z/ G C and becomes 
the identity when u is real. The formula for the new curve depends in general 
on two complex parameters: and the parameter c = c_|_/c_, assuming x = 
c+x"*" + C-X~ is expressed on fixed basis {x^ ■iX~) of solutions of the linear system 
at (k, v). 

2.3 Closure conditions 

When discussing the conditions for which the curve produced by a double 
Backlund transformation (36) is closed, we distinguish two cases: 

Case (A) is analogous to the situation discussed in §1.5 for a single Backlund 
transformation. For an initial curve 7 of length L, x(-s; ^) is taken to be an 
eigenfunction of the transfer matrix ^{kL, v) across k periods of 7. (However, 

when 7 is odd and k is odd, the appropriate transfer matrix is ^° q) ^(^-^)^) 
instead.) Then formula (36) produces a family of curves of period kL parametrized 
hj u eC. 

However, if the initial curve is a closed elastic rod of constant torsion, then one 
can show directly that, for any e C, the new curve has length L and is congruent 
to the original elastic rod. To see this, first note that the Riccati equation for z is 
now 

= z^cna;+ -(1 - ^^). (38) 

Since, when k — cnx, '^{kL; —ia/2) has eigenvectors given by (13) when o" e M, 
similar formulas hold for the eigenvectors of ^{kL; v): 



v+ = A 



(7 



2/x 

V V 2// 



v- = A 



(39) 



\ 2ii I 



where, now, a = 2iu for arbitrary u & C, and Ai and p are defined in terms of 
a as in §1.3. We again observe that, since the eigenvectors are independent of k, 
the Backlund transformation has period kL if and only if it has period L. We will 
consider the initial value problem for (38) corresponding to the other being 
similar; the initial value for z is then 

2a -4p 

-2(0) = o O) 
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which can be rewritten as 

sn a 

z(0) = ipsna — dna, where = 2pi'. (40) 

1 + cna 

Proposition 2.4. Let p G (0, 1) be the elliptic modulus, and let a lie in an open 
disk about the origin in the complex plane such that en (a) —1. The solution of 
the initial value problem (38), (40) is given by 

z+ = — (dn X + ipsa x) (dn {x — a) + ipsa {x — a)) . 

Moreover, Zj^{x) satisfies the identity 

cnx + ii' ( ] — cn{x — a). (41) 

Proof. For a, x real, z must be of unit modulus, and the substitution z = —e~^^ 
yields the initial value problem (14), (15) solved in §1.5. The solution is 

P+{x) = — [sin~"^(psna;) + sin~"^(psn (x — a))] 

and z+{x) — — e~*^+ has the above form and satisfies the required identity. The 
rest follows by analyticity in a and x. 

After a single iteration of the Backlund transformation, we have k^^^ — en (x—a) 
by the above identity. Let y = x — a. After the double Backlund transformation, 
we have 

where is a solution of 

-- = z.cny+-(l-.). 

The formulas in §2.2 show that, in performing the second Backlund transformation, 
we use the solution w+ = z+ ( — — ] . We also know from Prop. 2.4 (with 



X changed to y and v to P) that 

w = — (dn J/ + ipsny) (dn {y — d) + ipsa {y — a)) 

is a solution. In fact, it can easily be verified, for example by evaluated at a; = 0, 
that w and w+ are the same. Now it follows by the identity (41) that k^^^ = 
cn (a; — a — a), and the new curve is congruent to the original elastic rod. 

Case (B). When the transfer matrix is a multiple of the identity matrix, x(s, ^) can 
be taken to be an arbitrary complex linear combination c+x"'"(s, v) + c_x~(s, v) 
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of solutions to (26). To see when this is possible, let A(/cL, z^) be the trace of 
the transfer matrix; in Floquet theory, this is known as the Floquet discriminant 
associated to the linear system. 

Assume for the moment that k is even or 7 is even. Since det$(/cL, z/) = 1, 
the eigenvalues of the transfer matrix will be both 1 or both —1 exactly when 
= 4. Moreover, if 1/ is a root of multiplicity two for the function — 4, then 
the corresponding eigenvectors will be linearly independent. Taking into account 
the obvious changes for the case when k is odd and 7 is odd, we summarize the 
discussion in the following 

Proposition 2.5. Let u be a complex double root of the equation (a) A^(A;L, u) — 
4 = 0, or (b) /S?{kL, //) + 4 = i/7 is odd and k is odd. Let = c_|_x"'"(s; v) + 
C-X~{s\p), where x^,x^ oire the columns of the fundamental matrix of (26). 
Then, formula (36) produces a family of closed curves of period nL, parametrized 
hy uj = Cj^ / C- e C. 

We will loosely refer to the new closed curves as A;-fold covers of the original 
curve. 

We should note that, when the initial curve is an elastic rod of constant torsion, 
and the initial value for x is equal to the limit, as z/ ^ i/, of either Vj^ or V-, then 
the arguments in Case (a) apply to show that the new curve is congruent to the 
original elastic rod. 

In the next section we will exhibit concrete examples of new closed constant 
torsion curves obtained by double Backlund transformations of elastic rods of 
constant torsion. 

2.4 Exotic Curves from Elastic Rods 

To produce new curves, we need to solve the equations in Prop. (2.5) for the 
appropriate values of v. Suppose as before that the rod is of length L and this 
corresponds to the elliptic parameter x running from zero to 2nK. Then the 
formulas in §1.4 give 

(tr$(/cL))2 - 4 = ir{^{kLf) - 2 = -4 sin^ (/cnKA) 

when kn is even, and when kn is odd, 

(tr (; J) $(fcL))2 + 4 = tr(((j ^{kL)f)+2^ 4sm\knKA) 

(Note that while the formulas in §1.4 apply to the linear system (5), the fun- 
damental matrix for (26) differs only by conjugation by matrix A, and so the 
trace formulas give the same expression.) So, in either case we look for zeros of 
sm{knKA), where A is defined by (9) as a function of p and a, and a = 2iv as 
above. 

In practice, we look for values of cr in the first quadrant, because changing f 
to V does not change the double Backlund transformation, while changing v to 
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-2.5 2.5 5 7.5 



Figure 2.1. Knotted double Backlund transformation of the unknotted 
(1,3) elastic rod {p ^ .63093, n = 3, /c = 4, a ^ 1.2283 + 0.9688z). In 
the oblique views, the curve itself is the centerline of the vertical ribbon; 
the elastic rod is at left, the transformed curve at right. In the view 
from above, the curve itself is the inner edge of the radial ribbon. 

—u is equivalent, because of symmetries in the linear system, to just changing 
the sign of t^; = c_(_/c_. The formula (10) shows that sin^(-ftrA) has simple zeros 
at 0" = ±p' /p ± i — that is, where n = a/ (p~^ — cr^)^ + 4a'^/4 vanishes. (In fact 
these seem to be the only complex zeros of sm'^{KA); there are, however, an infinite 
number of real zeros, since A is asymptotically linear in a when a G M (cf. formula 
(10)).) We find, using Newton's method, simple zeros for sin(A:n-ftrA) that have 
real part close to p' /p. Once such a root has been found, we can generate a family 
of double Backlund transformations parametrized by a; G C; unless specified, the 
reader may assume the value used is w = 1. 
We observe the following phenomena: 

Number of roots increasing with k: Of course, all the zeros for sm[knKA) remain 
when we multiply k by an integer. However, zeros are also picked up from sm{iKA) 
where £ is a divisor of kn. For example, when the original curve is the closed 
unknotted elastic rod of constant torsion corresponding to p ~ 0.63093 and n = 3 
(appearing in Figure 2.1), there is one root for k = 2, two new roots for /c = 4, 
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and four new roots for k = 8. 

Obtaining knots from unknots: Performing a double Backlund transformation of 
the aforementioned unknot, using one of the roots for A; = 4, yields the knotted 
curve of constant torsion shown in Figure 2.1. One can verify that this knot has a 
minimum crossing number of 12, and is not a torus knot: its Alexander polynomial 
is 

A{t) = l-t + t^ -t^ + t^ -t^' + t'' -t^ + t^"", 

while the only torus knot that has an Alexander polynomial of this degree is a 
(2, ll)-knot, whose polynomial is different. So, this is a new knot type realizable 
by curves of constant torsion. 

Knotting is related to parametric resonance: When, for a fixed initial curve and 
fixed value of k, there are several roots available, choosing the a- value closest to 
p' /p + i produces the Backlund transform that has the most complicated shape 
(see below for examples). At a = p' /p + i, the Floquet multipliers coincide but 
the transfer matrix has only a single eigenvector up to multiple, and the general 
solution of the linear system is unbounded. 

Obtaining unknots from knots: When the original curve is the elastic rod of 
constant torsion which gives a (2, 5) torus knot, corresponding to p ~ 0.7845 and 
n = 5, three different roots are available for k — 2. Two of the resulting curves, 
shown in Figure 2.2, are unknotted. 

Self-intersections and dependence on uj: The shape of the transformed curve seems 
to depend more on the argument of u than its magnitude. For uj chosen to be real, 
we sometimes find that the transformed curve has self-intersections. For example, 
the third available root for the (2, 5) knot with k = 2 gives, using w = 1, the self- 
intersecting curve with 180-degree symmetry in Figure 2.3. The self-intersection, 
along with the symmetry, disappears when uj = e*"^/^. The resulting embedded 
curve is the twelve-crossing knot shown in Figure 2.3. The shape of the curve 
changes continuously as the argument of co varies, showing that a given knot type 
can have non-congruent realizations as a closed curve of a given constant torsion. 
Moreover, this curve has the same knot type as the curve in Figure 2.1, showing 
that this knot type can be realized by closed curves of different constant torsion! 
(Recall that the Backlund transformation preserves the value of the torsion, so 
these curves have the same torsion as the elastic rods we started with.) 

Finally, we note that these exotic curves, produced by double Backlund trans- 
formations of closed elastic rods of constant torsion, are planar- like 3-solitons — i.e., 
a linear combination of the purely binormal LIE vector fields Bq, Bi, B2, B^ van- 
ishes along them. This follows from our mechanical verification, for n = 0, 1, 2, of 
the conjecture in §1.7, since the calculation is purely algebraic and applies when 
the Backlund parameter C takes successive complex values u and u. 

2.5 Further Research 
In this section, we outline some directions for future research. 
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Figure 2.2. Unknotted double Biicklund transformations of the knot- 
ted (2,5) elastic rod, with elastic rod shown at above left {p ~ .7845, 
n = 5, /c = 2; a ^ 0.8982 + 0.8714z above right, a ^ 0.8821 + 0.6716z 
below) . 

Evolution under the LIE: It would be interesting to see how our exotic curves 
evolve under the Localized Induction Equation: 

Since kB is, up to a tangential component, a Killing field along the elastic rods, 
they evolve by a rigid motion. (In this case, the LIE preserves constant torsion; we 
do not expect this to happen in general.) In this context, Ercolani et al [E-F-M] 
used Backlund transformations of a given solution u{x, t) of the sine-Gordon equa- 
tions to produce solutions that appear as homoclinic orbits, as t goes from — oo 
to oo, in the isospectral set of u. For the LIE itself, the first author [C] has pro- 
duced homoclinic solutions as Backlund transformations of the planar circle (which 
translates under LIE); moreover, these solutions, which are multiple covers of the 
circle, exhibit self-intersections that persist under LIE. It is possible that some of 
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Figure 2.3. Resolving self-intersection for a double Backlund trans- 
formation of the (2,5) elastic rod [p ~ .7845, n = 5, /c = 2, a ~ 
0.9067 + 0.9697z). The upper curve has two self-intersections, one of 
which is visible in the closeup on the right. The lower curve differs by a 
change of initial values; at right is a view of it from above. 



the self-intersecting LIE solitons produced by our Backlund transformations may 
also exhibit this behaviour. 

The space of constant torsion knots: Our examples of exotic curves of constant 
torsion raise the question of which knot types can be realized with constant torsion. 
For example, while the trefoil knot is not realized by constant torsion elastic rods (it 
can, however, be realized by elastic rods of non-constant torsion [I]), is it possible 
that by a combination of knotting and unknotting, effected by successing double 
Backlund transformations of such elastic rods, a trefoil can be produced? More 
generally, why, in all of our examples, do all the crossings have the same sign? 

Understanding knotting: The examples given in the previous section show that, 
under the double Backlund transformation, the knot type of the curve may change 
dramatically. One approach to characterizing the change in knot type would be 
to find invariants which can be computed using the analytic representations we 
have for these curves. For example, the Mobius energy of the curve, which has 
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been shown to be related to the minimum crossing number [F-H-W], could be 
computed using the Gauss integral. It may be possible to compute the self-linking 
number of the new curve and the linking number of the new ciirvc and the old 
curve. (However, White's formula for the linking number is only valid when the 
ribbon stretching between the curves is embedded; we were able to use it in §1.6 
by taking the Backlund parameter C = to be small, but this not true for the 
roots V required for the double Backlund transformation.) It may also be possible, 
by generating and examining many more examples, to formulate conjectures con- 
cerning the change in knot type, in terms of constructions more familiar to knot 
theorists (eg. cable knots and satellite knots). 

We have shown how the Floquet spectrum of the linear system for an elastic 
rod contains information about both its knot type (related to the simple roots 
± i of the Floquet discriminant) and the knot types of its iterated Backlund 
transformations (encoded in some way by the complex double roots). Another 
approach, then, is to use Floquet analysis in combination with Backlund transfor- 
mations for general multi-phase solutions. Exact formulas for curves of constant 
torsion that generate n-phase solutions of the LIE or the sine-Gordon equation can 
be constructed in terms of Theta functions using standard methods of algebraic 
geometry [Kr, D] which, at the same time, provide a description of the associated 
Floquet spectrum. (Counting-type of arguments [P-T] can also be used to deduce 
a priori information about the Floquet eigenvalues.) 

Filling up the space of knots: One long-term goal is a complete classification of the 
knot types of n-phase solutions. In the spirit of Gromov's knot approximation by 
Legendrian curves [G] , it may be possible to combine density theorems for n-phase 
curves with a precise knowledge of their topology in order to approximate more 
general knotted curves. 
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